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ABSTRACT 



Given the considerable percentage of stars that are members of binaries or steUar multiples 
in the Solar neighborhood, it is expected that many of these binaries host planets, possibly even 
habitable ones. The discovery of a terrestrial planet in the a Centauri system supports this 
notion. Due to the potentially strong gravitational interaction that an Earth-like planet may 
experience in such systems, classical approaches to determining habitable zones, especially in 
close S-Type binary systems, can be rather inaccurate. Recent progress in this field, however, 
allows to identify regions around the star permitting permanent habitability. While the discovery 
of a Cen Bb has shown that terrestrial planets can be detected in solar-type binary stars using 
current observational facilities, it remains to be shown whether this is also the case for Earth 
analogues in habitable zones. We provide analytical expressions for the maximum and RMS 
values of radial velocity and astrometric signals, as well as transit probabilities of terrestrial 
planets in such systems, showing that the dynamical interaction of the second star with the 
planet may indeed facilitate the planets detection. As an example, we discuss the detectability of 
additional Earth-like planets in the averaged, extended, and permanent habitable zones around 
both stars of the a Centauri system. 



Subject headings: Planets and satellites: detection 
analytical 

1. Introduction 

The past decades have seen a great number 
of discov eries of planets aroun d stars other than 
our Sun ( Schneider et al. 201ll) . As some of these 
planets are of terrestrial nature, hopes to identi- 
fying Earth analogues have lead to considerable 
advances to wards the detection of possibly habit- 
able worlds ('Bot^IIMI l^rdeE^lliQl^ Even 
though quite frequent in the Solar n eighborhood 
( Kiseleva-Eggleton fc Eggleton 200 not many 
attempts have yet been made to specifically tar- 
get binary stars in this endeavor. Nonetheless, 
more than 60 planets have already been found 
in and around such systems ( Haghighipour 2010t 
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20121 ). Although several P-Type (circumbinary) 
planets orbiting b oth stars of a close binary have 
also been found (iDovle et al] 1201 11 I Welsh etldl 
20121 : lOrosz et al Il2012allbl). most planets are in 
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the so-called S-Tvpe (|Rabl fc Dvorak 1 19881 ) con- 
figuration where the planet orbits only one of the 
binary's stars. A prominent example of an S-Type 
system is a Centauri AB which hosts a terrestrial 
planet around the fainter bin ary component, a 
Cen B (jPumusaue et al. II2012I ). 

The reason for the general reluctance to in- 
clude binary systems in the search for terres- 
trial, habitable planets lies in the assumption that 
the additional interactions with a massive com- 
panion will make planets harder to find. That 
is primarily because the gravitational interaction 
between the second star and a planet may al- 



ter the planet's orbit significantly and compli- 
cate the task of interpreting the planetary sig- 
nal. One aim of this work is, therefore, to show 
that changes in the planet's orbit can actually en- 
hance its detectability (see section [S|). Of course, 
the orbit of a binary as well as its stellar pa- 
rameters have to be well determined in order to 
be able to identify signals from additional terres- 
trial planets. Sensing the need for a better un- 
derstanding of binary star systems, efforts have 
been intensified to improve physical as well as or- 
bital data for nearby binaries 
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Understanding the complex interactions be- 
tween a stellar binary and a planet is essential if a 
system's potential habitability is to be evaluated. 
For instance, one of the main assu mptions of clas- 



sical h abitability, as introduced by iKasting et al 



(|1993I ). is that the planet moves around its host 
star on a circular orbit. This may not be a valid as- 
sumption for planets in a binary star system where 
the gravitational perturbation of the secondary 
can excite the eccentricity of the pl a net's orbit 
(iMarchaj 'l990; 'Ceorgakarakosj l2002t lEggl et al ' 



2012t ) lEggl ct al. (2012) found that except for 
S-Type systems where the secondary star is much 
more luminous than the planet's host star, varia- 
tions in planetary orbit around the planet-hosting 
star are the main cause f or cha nges in insolation. 
Even though lEggl et al.l (l2012h gave an analytic 
recipe for calculating the boundaries of the HZs 
in S-Type binaries, it remains to be seen whether 
an Earth-like planet in the HZ of a system with 
two Sun-like stars will in fact be detectable. 

In order to answer this question, we consider 
three techniques, namely, radial velocity (RV), as- 
trometry (AM) and transit photometry (TP), and 
discuss whether the current observational facili- 
ties are capable of detecting habitable planets in 
such systems. We provide analytical formula for 
estimating the strength of RV and AM signals for 
habitable. Earth-like planets, and show that the 
planet-binary interaction can enhance the chances 
for the detection of these objects. 

The rest of this article will be structured as fol- 
lows. In sections [2] and |3] analytic estimates of 
the maximum and root mean square (RMS) of 
the strength of an RV and an AM signal that 



an Earth-like planet produces in an S-Type bi- 
nary configuration will be derived. Section |4] will 
deal with consequences of such a setup for TP. 
We will then briefly recall the differe nt types of 



HZs fo r S-Type binaries established in lEggl et al 



(|2012l ) , and use their methodology to identify sim- 
ilar habitable regions in the a Centauri system 
(section [5]) . This system has been chosen because 
first, it has inspired many studies on the possibility 
of the formation and detection of habitable plan- 
ets aroun d its stel lar components (iForga nl l2012t 
Guedes et al. 20 0l iThebault et al.ll2009[ ) and sec- 



ond, Dumusque et al. I ( 2012[ ) have alreadv discov- 
ered an Earth-sized planet in a short-period orbit 
around its secondary star. Therefore, we will com- 
pare our RV estimates to the actual signal of a Cen 
Bb, and study its influence on an additional ter- 
restrial planet presumed in a Cen B's HZ. Finally, 
in section m the projected RV, AM and TP traces 
that terrestrial planets will leave in the HZ of the 
a Centauri system are analyzed, and the results 
are discussed within the context of the sensitivity 
of the current observational facilities. 

2. Radial Velocity 

To estimate the RV signal that an Earth-like 
planet produces in an S-Type binary system, 
we will build upon the formalism presented by 
iBeauge et al. ( 20071 ). We assume that the non- 
planetary contributions to the host star's RV sig- 
nal (such as the RV variation caused by the mo- 
tion of the binary around its center of gravity) are 
known and have been subtracted, leaving behind 
only the residual signal due to the planet. The 
motion of the planet around its host star then 
constitutes a perturbed two body problem, where 
the gravitational influence of the secondary star 
is still playing a role and is mirrored in the forced 
variations of the planet's orbit. 

In practice, the extraction of the planetary sig- 
nal is all but a trivial task. Even after subtrac- 
tion of the binary's barycentric and proper mo- 
tion, the residual will contain contributions from 
the binary's orbital uncertainties as well as from 
non-gravitational sources which could be orders 
of magnitude larger than the star's reflex signal, 
such as the Rossiter-McL aughlin effect in t ransit- 
ing systems, for example (jOhta et al.l 120051 ). The 
discovery of a Cen Bb showed, however, that a 
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substantial reduction of non-planetary RV inter- 
ference is possible if the respective binary star has 
been studied in sufficient detail. 

The amplitude of the planet induced RV signal 
of the host star, Vr, is given by 

Vr = K [cos(/ + uj) + ecoso;] , (1) 
where K is equal to 



K 



/i (k np) 



1/3 . 



(2) 



In equation ([2]), ^, — mi/ {mo+mi) is the planet to 
star mass-ratio with mi and mo being the masses 
of the planet and host-star, respectively. The 
planet's mean motion, Up — 27r/Pp, is given by 
Up = \J njo? with K — t/(mo -I- mi), and Pp and 
Q being the planet's orbital period and gravita- 
tional constant. The quantities a,e,i,f and uj in 
equations (1) and (2) denote the planet's semima- 
jor axis, eccentricity, orbital inclination relative to 
the plane of the sky, true-anomaly, and argument 
of periastron, respectively. 

Our goal in this section is to identify the range 
of the possible peak amplitudes that a terrestrial 
planet in an S-Type binary configuration can pro- 
duce. We note that the gravitational infiuence 
of the second star causes the planet's orbital el- 
ements to vary, thus inducing additional time de- 
pendent changes in the radial velocity signal Vr 
( Lee fc Peald 120031 ). While we know from secular 
perturbation theory that a does not change sig- 
nificantly with time for hierarchic al systems such 
as the one under c onsideration ( Marchail 1990t 
ICeorgakarakos 2003), w becomes a function of 
time. We assume coplanar orbits of the planet and 
the binary star which results in the planet's incli- 
nation to the plane of the sky (i) to remain con- 
stant. In contrast, the planetary eccentricity will 
vary between zero and a maximum e"*"'^, where 
the latter value can be expressed as a function 
of the syste m's masses and the binary's orbital 
parameters ( Eggl et al. 20121) . This is important, 
because the reflex RV signal (Vr) of a star can be 
increased significantly by planetary orbital eccen- 
tricities (Fig. [Ij. Using equation ([1]) we identify 
the global maximum of at / = w = 0, when 



= V 



where 



nil smi 



^a(mo + mi) 



(3) 



(4) 



Equation (j3|) presents a fully analytic estimate of 
the expected maximum RV signal that a terres- 
trial planet produces in an S-Type binary config- 
uratior0. 

As an example for the influence of a double star 
on a planetary RV signal, the induced variations 
in the RV of the planet's host-star are presented in 
figure [H The host-star is a constituent of a Solar- 
type binary with a semimajor axis of 20 AU and 
an orbital eccentricity of 0.5. Changes in the am- 
plitude of Vr are due to variations in the planet's 
eccentricity. 

Since we do not know the state of the planet's 
orbital eccentricity at the time of observation, we 
consider a range for the maximum possible ampli- 
tudes of its radial velocity. 



T/circ ^ I ^ J /max 

Vr < Vr\f=Lj=Q < Vr 



(5) 



Although the range of the amplitude of the host 
star's RV signal, as given by equation ([5|), can be 
used to identify the "best case" detectability lim- 
its, the maximum values of the RV signal due to 
the planet will be "snapshots" that are reached 
only during brief moments. As a result, their val- 
ues for assessing the precision needed to trace fin- 
gerprints of an Earth-like planet are rather lim- 
ited. In such cases, expressions for the RMS of 
the astrometric signal are preferable. 

Since RMS values are by convention time- 
averaged, we substitute / by the mean anomaly 
M — Hpt in all corresponding functions using the 
equation of the center expansion up to the sixth 
order in planetary eccentricities (see appendix E| 
and average over M and uj. The vastly different 
rates of change of these quantities (M ^ uj) make 



^Larger signals are possible, if the terrestrial planet 
has a considerable initial eccentricity after its forma- 
tion and migration phase. Yet, due to the eccentric- 
ity dampening in protoplanetary discs this seems unlikely 
iPaardekooper fc Leinhardtll2010l ). 



3 



it possible to consider w to remain constant dur- 
ing one cycle of M , so that independent averaging 
can be performed. In order to eliminate short term 
variations in the RV signal, we first average over 
M . Averaging over w as well might be desirable 
if for example the initial state of uj is unknown, 
or if observations stretch beyond secular evolution 
timescales of the planets argument of pericenter. 
We, therefore, define two different types of RMS 
evaluations for a square-integrable function F] 



planetary orbital periods are available in the data 
set. In the case of a Cen Bb, there are order-of- 
magnitude differences in the rates of change of the 
mean anomaly (M) and the argument of pericen- 
ter (w). It would, therefore, make more sense to 
assume w to be constant and add it as a variable 
in the fitting process. If stronger perturbations or 
additional forces act on the planet, the periods can 
be considerably shorter, so that the fully averaged 
equations might come in handy. 



{{F))m = {F^f^^ 



and 



M.i 



2ti 



27r 



2ir 



27r 



F^{M)dM 



1/2 



(6) 



F^{M,uj)dMduj 



1/2 



(7) 



3. Astrometry 

In order to derive the maximum and RMS val- 
ues for an astrometric signal, we will u se the 
framework presented in IPourbaixl (l2002h . We 
again assume that the non-planetary contributions 
have been subtracted from the combined signal of 
the host star and planet. The projected motion of 
the planet on the astrometric plane is then given 

by 



Using equations ^ and (O , the RMS values of Vr 
are then given by 



{(Vr)) 



with 



M 



{{Vr))M. 
1 - 



+ 0((e^)L) 



cos(2aj) 



1/2 



xe — cos E — 
He = B( cos E — e) 



F\/l - e2 sinE, 



GVi - e2 sinE, 



(10) 



where xe and i/e are the Cartesian coordinates of 
(8) the projected orbit, e is the planet's orbital eccen- 
tricity, E is the eccentric anomaly, and A, B, F 
and G are the modified Thiele-Innes constants 
given by 



{{Vr))M.u 



mil sin i I J_yc^rc ^g) 



yj2a{mQ + nil) 



Here we have considered {a)M = a since a ~ 
(|Marchallll990tlGeorgakarakosll2003l) . Also, 



{e^)M cos(2a;) duo — 0, 



as indicated in appendix |B1 It is noteworthy that 
the averaging over lu causes the RMS value of the 
RV signal to become independent of e so that its 
difference with the peak signal in the circular case 
fycirc\^ becomes a mere factor of 1/72. Thanks 
to their intricate relation to power-spectra, RMS 
values can also be valuable for orbit-fitting. The 
choice of singly or doubly averaged RMS rela- 
tions for this purpose will depend on how many 



A= — (coswcosfJ — sinw sini7 cosi), 

B= — (cos w sin f2 + sin w cos SI cos i), 

a (11) 
F= — — (sinwcosfi + cosw sinfJ cosi), 

G = — — (sin UJ sin SI — cos uj cos SI cos i) . 

In these equations, d is the distance between the 
observer and the observed system in units of the 
planetary semimajor axis a. We can rewrite equa- 
tions (|10l) in terms of the true anomaly / as. 



Xf = — r cos / H r smj, 

a a 

^ f^^ ■ f 
Uf = — r cos / H r smj. 



(12) 



a 



a 
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In these equations, r — a(l — e^)/(l + ecos/) 
represents the planet's radial distance to its host 
star. Because the motion of the planet itself can- 
not be traced, we translate these equations into 
the apparent motion of the host star by the appli- 
cation of Newton's third law. That is, 



— X — fix f , 



(13) 



Here, X and Y are the projected coordinates of 
the center of mass of the planet-star system, and 
IJ, denotes the planet-star mass-ratio as defined for 
equation (2). 

Assuming without the loss of generality that 
the barycenter of the star-planet system coincides 
with the origin of the associated coordinate sys- 
tem, the distance of the projected stellar orbit to 
the coordinate center will be equal to 



2 2 



l?a^ (1 



e2)2[l-sin^sin2(/ + w)] (14) 



(1 + ecos fY 



The right-hand side of equation (|14p is inde- 
pendent of and has a global maximum at 
f — TT^Lj — Q when e = ^-max ^ This translates into 
a maximum astrometric amplitude given by 



where 



(1 



0, 



(15) 



(16) 



The planetary maximum AM signal will again lie 
between and p™"^. A remarkable feature of 
piriax pcirc their independence of the sys- 
tem's inclination i. This is visualized in figure [21 
The same figure also shows the time evolution of 
the AM signal due to an Earth-like planet orbiting 
a Cen B at a distance of 1 AU. 

The astrometric RMS values are given by 



M 



3(e2 



>M 



(5 cos[2a;] — 3) I sin' 



1/2 



and 



M,i 



9 



^ + )cos(2i) 



1/2 



(18) 



Details regarding the derivation of equations ([T7|) 
and ([T8| can be found in appendix [BJ In contrast 
to the doubly averaged equations for the RMS of 
an RV signal, equation ^TE\\ shows a dependence on 
the binary's eccentricity. In cases where the plan- 
etary inclination i coincides with the inclination of 
the binary its elf, analytic expressions fo r {e'^)M,Lj 
are available (|Georgakarakosll2003l l2005lf l. 



4. Transit Photometry 

In transit photometry, signal strength is equiv- 
alent to the relative depth of the dint the planet 
produces in the stellar light-curve during its tran- 
sit. Assuming that the star-planet configura- 
tion allows for occultations, and excluding graz- 
ing transits, the fractional depth of the photomet- 
ric transit (TD) produced by an Earth-like planet 
is simply given by the proportion of the lumi- 
nous area of the disk of the star that is covered 
by the planet as the planet moves between the 
observer and the star. Ignoring limb darkening, 
that means, TD ~ Rp/Rl where Rp is the radius 
of the planet and Ri, is the stellar radius. The 
overall probability to obser ve a transit is given by 
( Borucki fc Summers 11984 : 



R. 

Pt ^ — 



(19) 



In equation p9|) . rx is the radial distance of 
the planet to the star during the transit. For 
an eccentric planetary motion, the planet-star 
distance during transit can be expre s sed a s 



a(l-e2)/(l + 



e cos (jj 



(iFord et all l2008[ ). 



where ui denotes the argument of pericenter mea- 
sured from the line of sig In analogy to sec- 
tions [5] and [31 the maximum and averaged transit 



^The analytic expressions given in these articles are also 
averaged over initial phases, i.e. different relative starting 
(17) positions of the planet and the binary stars. 

•^Note that this is different from the conventions used for RV 
and AM measurements. 
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probability for a planet perturbed by the sec- 
ondary star in a planar configuration can be caf 
culated by substituting for tt in equation ()19| 
and averaging over uj. This will result in 



and 



Pt 



(J>t)lc 



ail 



a 1 



(20) 



(21) 



{e'^)M,uj) ' 

Equations (I^H)) and (pij) indicate that the increase 
in the eccentricity of the planet due to the pertur- 
bation of the secondary increases the probability 
of transit. In deriving these equations, we have 
ignored the occultation of the planet by the sec- 
ond star. However, depending on the period ra- 
tio between the secondary and the planet, such 
conjunctions are either scarce or short-lived. Con- 
sequently, their contribution to the probability of 
witnessing a planetary transit is negligibly small. 

5. Application to the a Centauri system 

In this section we will show that the previously 
derived analytic expressions produce results that 
are in good agreement with the current observa- 
tions of a Cen Bb. We will also present numerical 
evidence that the presence - or absence - of an 
additional terrestrial-planet in the HZ of a Cen 
B cannot be derived easily from the orbit evolu- 
tion of a Cen Bb. Consequently, we argue that 
an independent detection of additional terrestrial 
companions might be difficult, but more promis- 
ing. For this purpose, we will determine the HZ 
of a Cen B, as weU as the RV, AM and TP sig- 
natures of an Earth-like planet orbiting in the HZ 
of a Cen B. Since there is no a priori reason why 
the brighter component of a Centauri could not 
be hosting a terrestrial planet as well, we perform 
a similar analysis for a Cen A. We will also study 
the behavior of equations p1l2T|) for a broad range 
of binary eccentricities. 

5.1. a Centauri's terrestrial planet 

The planet discovered around a Cen B offers a 
perfect opportunity to compare the RV amplitude 
predictions derived in section [2] with actual mea- 
surements. The planet's known orbital parame- 
ters are given in Table [TJ In Table [2] we present 



the analytic estimates of section [2] applied to the 
a Centauri ABb system. Assuming the system to 
be coplanar {i ~ 79.2°), the predicted RV ampli- 
tude for circular planetary motion (V^"''') is very 
close to the observed RV amplitude. This is not 
surprising, since the planetary parameters were 
derived from an RV signal using the same method- 
ology in reverse. While still well within measure- 
ment uncertainties, the deviation of the maximum 
RV amplitude {V^""^) from the observed value is 
larger than that of V^'^'^'^. On the one hand, this 
might indicate that the planet is currently in an 
orbital evolution phase where its eccentricity is al- 
most zero. On the other hand, the planet may 
be too close to its host star for our model to pre- 
dict VT""^ correctly. In fact, we show in section 
I5.3l that the latter explanation is more likely, since 
the influence of general relativity (GR) cannot be 
neglected in this case. Estimates based on Newto- 
nian physics exaggerate the actual eccentricity of 
a Cen Bb. Its orbit remains practically circular 
despite the interaction with the binary star (see 
section [531 for a detailed discussion). This justi- 
fies the assumption of a c i rcular planetary orbit 
made bv lDumusque et al~ ( 20121 ). 

Since we are especially interested in additional 
habitable planets, however, it is worthwhile to ask 
whether predictions on the orbital evolution of a 
Cen Bb can be used to exclude the presence of 
other gravitationally active bodies in the system. 
In other words: Could an Earth-like planet still 
orbit in the HZ of a Cen B or would the accom- 
panying distortions of the orbit of a Cen Bb be 
significant enough to detect them immediately? 
Before we try to answer these questions, we need 
to briefly recall some important aspects regarding 
HZs in binary star systems. 

5.2. Classification of HZs 

Combining the c lassical definition of a HZ 
( Kasting et al. 1993[ ) with the dynamical prop- 
erties of a plan et-hosting double star system, 
Eggl et al.l (|2012l ) have shown that one can dis- 
tinguish three types of HZ in an S-type binary 
system: 

The Permanently Habitable Zone (PHZ) 

where a planet always stays within the 
insolation liin i ts (Sj , S g) as defined by 
Kasting et al. ( 19931) and lUnderwood et al 
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(j2003l ). In other words, despite the changes 
in its orbit, the planet never leaves the clas- 
sical HZ. The total insolation the planet 
receives will vary between the inner (Sj) 
and outer (So) effective radiation limits as 
Si > Stot > So where, for a given stellar 
spectral type, Sj and So are in units of 
Solar constant (1360 [W/m^]). 

The Extended Habitable Zone (EHZ) where, 
in contrast to the PHZ, parts of the plan- 
etary orbit may lie outside the HZ due to 
the planet's high eccentricity, for instance. 
Yet, the binary-planet configuration is still 
considered to be habitable when most of the 
planet's orbit remains inside the boundaries 



of the HZ. In this case, {St, 



a < Si 



and {Stot)t~'^ ^ So where {Stot)t denotes 
the time-averaged effective insolation from 
both stars and is the effective insolation 



variance. 



The Averaged Habi table Zone (AHZ). Foll owing 
the argument of IWilliams fc PollardI (|2002h 
that planetary eccentricities up to e < 0.7 
may not be prohibitive for habitability as 
long as the atmosphere can act as a buffer, 
the AHZ is defined as encompassing all con- 
figurations which support the planet's time- 
averaged effective insolation to be within 
the limits of the classical HZ. Therefore, 
Si > {Stot)f > So- 
Analytic expressions for the maximum insolation, 
the average insolation {{Stot)), and insolation vari- 
ance that a planet e ncounters in a bin ary system 
have been derived in lEggl et al.l (|2012h . We refer 
the reader to that article for more details. 

Figures [3] and m show the application of the pro- 
posed habitability classification scheme to the a 
Centauri system. In these figures, blue denotes 
PHZs, green shows EHZs, and yellow corresponds 
to AHZs. The red areas in figures [3] and |4] are 
uninhabitable, and purple stands for dynamically 
unstable regions. Table[T]shows the physical prop- 
erties of the syst e m. W e used the formulae by 
Underwood et al ] (|2003[ ) to calculate Si and 5*0 
for the given effective temperatures of a Cen A 
and B. In general, these formulae allow for ex- 
te nding the analytic estimates for HZs, as given 
bv lEggl et ai] (|2012l ). to main sequence stars with 



different spectral types. Runaway greenhouse and 
maximum greenhouse insolation limits were used 
to determine the inner and outer boundaries of 
HZ, respectively. 

As shown in figures [3] and HI the locations of 
the HZs and the detectability of habitable planets 
in those regions depend strongly on the eccentric- 
ity of the binary (ef,). The actual eccentricity of 
the a Centauri system is denoted by a horizontal 
line at Cf, = 0.5179. The values for the borders 
of the different HZs using a Centauri's actual ec- 
centricity are listed in Table [H As shown here, 
both stars permit dynamical stability for habit- 
able. Earth-like planets. Due to the difference in 
stellar luminosities, the HZs around a Cen A are 
larger and farther away from the host star com- 
pared to a Cen B. Since the binary's mass-ratio is 
close to 0.45, the gravitational influence of a Cen 
B is more pronounced on the PHZ of a Cen A. 
This is a consequence of the larger injected plan- 
etary eccentricities (cp) as can be seen from the 
top row of figure [5j The relatively larger gravi- 
tational influence of a Cen B onto the HZ of a 
Cen A is also mirrored in the fact that the region 
of dynamical instability (purple) reaches towards 
lower binary eccentricities. The change in the 
range and conflguration of HZs with the change 
in planetary semimajor axis and eccentricity of 
the binary is pronounced. A clear shrinking trend 
for PHZ and EHZ can be observed for high val- 
ues of_thfi_bina£;Yls eccentricity. While as shown 
bv lEggl et ai] (|2012[ ). the AHZ in general expands 
slightly when the eccentricity is enhanced, figures 
[3] and |4] show that in the a Centauri system, this 
HZ depends only weakly on ei, making it the clos- 
est approximation to t he classical HZ as defined 
bv lKasting et al. ( 1993 ). Comparing these results 
with the existing stud i es on the H Zs for a Cen B 
such as Gucdcs et al.l (|2008[ ) and iForganl (|2012l ). 
one can see that the values of the inner bound- 
aries of the HZs around a Cen B as given in fig- 
ures PandB l coinci de well with the previous stud- 



ies. 



Forgan (l2012l) even found a similar shrink- 



ing trend with higher planetary eccentricity. Yet, 



Forganl (j2012l ) did not take the actual coupling 



between the planet's eccentricity and the binary's 
orbit into account. The limits for the outer bound- 
aries of HZ in our model are different from the ones 



in iForganI (|2012[ ) since different climatic assump- 
tions were made. In this work we used insolation 
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limits for atmospheric collap se assuming a maxi- 
mum greenhouse atmosphere ([Kasting et al 
whereas IPorgan ( 2012.) focused on emergence from 
snowball states. 

5.3. Additional terrestrial planets in a 
Centauri's HZs 



While the classification of habitable zones pre- 
sented in the previous section is globally applicable 
to binary star systems , the analytic est imates to 
calculate their extent ( Eggl et "al1l2012[ ) are only 
strictly valid for three body systems, e.g. the bi- 
nary star and a planet. Additional perturbers will 
influence the shape and size of the HZs. It is thus 
necessary to investigate which eff'ect the already 
discovered planet around a Cen B would have on 
an additional terrestrial planet in a Cen B's HZ. 

If the mutual perturbations were large, the HZ 
boundaries given in Table [3] would have to be 
adapted, but a Cen Bb's orbital evolution could 
also contain clues on the presence - or absence - 
of an additional planet. Should the interaction 
between the inner planet and an additional ter- 
restrial body in the HZ be small, then the HZ 
boundaries would hold. However, a detection of 
the habitable planet via its influence on a Cen 
Bb's orbit would become difficult. 

In figure |6l results of numerical investigations 
on the coupled orbital evolution of an additional 
terrestrial planet and a Cen Bb are presented. 
The top row of figure |6] shows the eccentricity evo- 
lution of a Cen Bb altered by an additional Earth- 
like planet at the inner (red curve) and outer (or- 
ange curve) edge of a Cen B's AHZ . The corre- 
sponding reference curve (blue) represents a Cen 
Bb's eccentricity influenced only by the binary a 
Cen AB. The top left panel of figure 3 shows the 
results in Newtonian three (3BP) and four (4BP) 
body problems. The top right panel depicts simi- 
lar analysis with general relativity (GR) included. 
The difference between the two approaches is quite 
pronounced, as GR clearly prevents the secular 
rise in Ce n Bb's eccentricity predicted in the classi- 



cal se tup (jBlaes et al.ll2002 



ty predictea m tne classi- 
3; Fabrvckv fc Tremaind 



20071 ). Thus, the orbit of a Cen Bb stays circular, 
even when tidal forces are neglected. The varia- 
tions in semimajor axis (Aa) for a Cen Bb are not 
shown, because they remain below 10~^ AU for all 
cases. 



A possible method to search for additional com- 
panions is to measure variations in a Cen Bb's or- 
bital period. Yet, the small Aa values make this 

1 /3 

approach diflScult, since AFp oc Fp'"Aa. Disen- 
tangling the effects of GR and perturbations due 
to other habitable planets on a Cen Bb's period 
would require precisions several orders of magni- 
tude greater than currently available. The top 
right panel in figurelHlshows that the perturbations 
an additional planet at the inner edge of a Cen B's 
AHZ causes in a Cen Bb's eccentricity (red) are, 
in principle, distinguishable from the nominal sig- 
nal (blue). Unfortunately, it is also clear from this 
graph that neither the required precision nor the 
observational timescales necessary to identify the 
presence of an additional Earth-sized companion 
via observations of a Cen Bb's eccentricity seem 
obtainable in the near future. For habitable plan- 
ets at the outer edge of a Cen B's AHZ the chances 
for indirect detection seem even worse, as their in- 
fiuence on a Cen Bb's orbit is negligible (orange 
vs. blue curves). 

In order to confirm that the interaction between 
a Cen Bb and Earth-like planets in the HZ is 
small, as well as to further study the infiuence of 
the GR on the dynamics of the system, we exam- 
ined the orbital evolution of a fictitious habitable 
planet in that region. The results are shown in 
the bottom row of figure [51 The left panel depicts 
the eccentricity evolution of additional terrestrial 
planets positioned at the inner (red) and outer 
(green) edges of a Cen B's AHZ. The secular vari- 
ations in the eccentricity (bottom left panel) and 
semimajor axis (bottom right panel) of the hab- 
itable planet were computed numerically, taking 
the infiuence of the binary a Cen AB, the planet 
a Cen Bb, as well as GR into account. When 
comparing the analytic estimates of e™"^ with the 
evolution of the habitable planet's eccentricity in 
the full system, it is evident that neither GR nor 
a Cen Bb alter the results for planets in a Cen B's 
HZ significantly. Also, the deviation in the hab- 
itable planet's semimajor axis due to GR and a. 
Cen Bb (Acp) remains below 0.1% and 0.5% for 
planets at the inner and outer edge of a Cen B's 
AHZ, respectively. 

We conclude that the interaction between ad- 
ditional terrestrial planets in a Cen B's HZ and a 
Cen Bb is indeed small. Thus, our estimates for 
the HZs of the a Centauri system remain valid. 
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The existence of additional terrestrial planets on 
the other hand cannot be determined easily from 
observing the orbital evolution of a Cen Bb. 

The presented results are, strictly speaking, 
only valid for a coplanar configuration, i.e. the 
binary and both planets are in the same orbital 
plane. Mutually inclined configurations can ex- 
hibit much more involved d ynamics such a s Koz ai 
resonant behavior (see e.g. Correia et al.l ( 201ll )). 
A detailed study of such effects lies beyond the 
scope of this work. Nevertheless, the arguments 
presented in this section suggest that the search 
for an additional coplanar planet in the HZ around 
a Cen B will most likely have to be performed 
without relying on observations of a Cen Bb. We 
will, therefore, investigate whether habitable plan- 
ets can actually be detected independently in Sun- 
like binary star configurations using current obser- 
vational facilities. 

5.4. Detectability Through Radial Veloc- 
ity and Astrometry 

We apply our methodology, as derived in sec- 
tions [2] and 131 to a fictitious terrestrial planet in 
the HZ of binary systems similar to a Centauri 
AB but with a broadened range of binary eccen- 
tricities. In addition to the habitability maps dis- 
cussed in section 15.21 figures [3] and |4] show the 
results regarding peak and RMS strength of the 
RV and astrometric signals. Here, the aim is to 
illustrate how the different types of HZs presented 
in section [521 as well as the maximum and RMS 
signal strengths defined in section [2l vary with the 
binary's eccentricity (et,) and planetary semimajor 
axis (Up). The left column of figure |31 shows maxi- 
mum (top) and RMS {bottom) values of the signal 
strengths for the more massive binary component, 
in this case similar to a Cen A. Results for the less 
massive component akin to a Cen B are shown in 
the right column. 

The dashed vertical lines in the top rows of fig- 
ures |3] and |4] represent the sections of the param- 
eter space with similar V^^'^'^ and p'^^^'^ values, re- 
spectively. Since V^^^'^ and p'^*'"'^ are independent 
of the planetary (and consequently the binary's) 
eccentricity, the different values of these quantities 
vary linearly with the planet's semimajor axis. In 
contrast, T^™"^ and p^"-^ ^ represented by the solid 
contour lines, depend on the maximum eccentric- 
ity of the planet (Cp"^^) and therefore change with 



the binary's eccentricity {eh). Since for circular 
binary configurations only small eccentricities are 
induced into the planet's orbit, V^"''^ and V^^^'^ 
almost coincide. The same holds true for p™°-^ 
and in this case. Yet, V™""^ and p'"'*^ grow 
with the binary's eccentricity. The corresponding 
contour lines indicate that for high binary eccen- 
tricities even small planetary semimajor axes can 
produce similar AM peak signal strength. Simi- 
larly, planets with larger distances to their host 
stars can still cause similar RV amplitudes if the 
binary's eccentricity is sufficiently large. If a fixed 
detection limit is set, e.g. Vr = 9.5 m/s, planets 
with semimajor axes up to 1.5 AU could still be 
found around stars similar to a Cen A, assum- 
ing a binary eccentricity of Cf, — 0.7. To pro- 
duce a similarly high RV amplitude, a circular 
planet has to orbit its host star at roughly 0.8 AU 
(Fig. [3]). In other words, high binary eccentricities 
lead to excited planetary eccentricities which in 
turn increase the peak signal strengths suggesting 
that binary-planet interactions can actually im- 
prove the chances for detecting terrestrial planets. 
Naturally, if the planet's eccentricity happens to 
be close to zero at the time of observation, this 
advantage is nullified. 

The bottom row of Fig. [3] shows the same 
setup with RMS signal strengths {{Vr))M,Lj and 
{{p))m,u,, respectively. While {{Vr))M,u, is inde- 
pendent of the binary's eccentricity, it is evident 
from equation ([T7| that {{p))m, uj depends weakly 
on Bb since (e^)^/^ ~ 0.1 for the cases considered 
and therefore (e^) <C 1 (see Fig. [5l bottom). The 
slight curvature of the contour lines representing 
the RMS signal in Fig. [H indicates this behavior. 
A summery of RV and AM signal strengths for an 
Earth-like planet at the boundaries of a Centauri's 
HZ is presented in Table [3l 

We illustrated in this section that the dynam- 
ical interactions between a terrestrial planet and 
the secondary star can produce large peak ampli- 
tudes which may enhance the detectability of the 
planet with the RV and AM methods consider- 
ably. The RMS values of the planet's AM and 
RV signals, on the other hand, remain almost un- 
affected by the gravitational influence of the sec- 
ondary star. 
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5.5. Transit Photometry 

To assess the detectabihty of a terrestrial planet 
in the HZ of a Centauri AB (and similar binaries) 
through transit photometry, we calculated the 
relative transit depths that an Earth-like planet 
would produce during its transit. If such a system 
hosted a transiting terrestrial planet, TD values 
would range around 55 ppm for a Cen A, and 
115 ppm for a Cen B. Such transit depths are de- 
tectable by NASA's Kepler telescope for instance 
- stellar and instrumental sources included - as the 
space craft's media n noise level amounts to ~ 29 
ppm (|Gilliland et a l. 2011). Therefore, Earth-like 
planets could in theory be found around a Cen- 
tauri stars. However, Kepler was not designed to 
observe stars with apparent magnitudes between 
and 3 such as those of a Centauri. The TESS 
mission (Transiting Exoplanets Survey Satellite), 
for instance, will aim for TP of brighter stars 
( Ricker et al.ll2010l ). Nevertheless, the example 
of Kepler suggests that the detection of transit- 
ing habitable planets in S-type systems would 
be possible using current technology. In fact, 
very much similar to the cases discussed in the 
previous sections, the eccentricity excitation that 
an Earth-like planet experiences in a binary star 
system may enhance its possibili ty of detection 
via transit photometry (Kane fc von Braun II2008L 
Kane. Horner fc von Braun I2OI2I : Borkovits et al 



20031 also see Fig. [7|). Assuming a Centauri was 
a transiting systenQ a comparison of the transit 
probabilities of actual planetary orbits to circu- 
lar orbits shows that an 18% increase in pT val- 
ues seems possible for terrestrial planets at the 
outer edge of a Cen A's AHZ (Fig. [T]). Given 
the right orbital configuration, it may be more 
likely to identify a transiting habitable terrestrial 
planet around a stellar component of a binary 
than around a single star assuming similar initial 
planetary eccentricities. 

The increase in transit probability for planets 
in double star systems is less dramatic when the 
equations are averaged over all possible configu- 
rations of the argument of pericenter as in equa- 
tion ([21]) . Averaged transit probabilities are rep- 
resented by the straight lines in figure [T] As 
{Pt) / {PT\e=o} > 1, for terrestrial planets' orbits 
with e > 0, the chance for transit is in general 



higher for Earth-like planets in binary stars than 
for terrestrial planets in circular orbits around sin- 
gle stars. 

6. Discussion 

Comparing the quantitative estimates of RV, 
AM and TP signals, TP seems to be the best 
choice for finding Earth-like planets in the HZs 
of a coplanar S-Type binary configuration with 
Sun-like components. Even for a system as 
near as a Centauri, AM peak signals only mea- 
sure /ias. Unfortunately, neither ESO's VLBI 
with PRIMA, nor ESA's GAIA mission will be 
able to deliver such prec ision in the near future 
(|Quirrenbach et al.ll201l[) . GAIA's aim to provide 
/ias astrometry will most li kely not be achieved un - 



til the end of the mission ( HestrofFer et al. 2010() 



Also, from an astrometric point of view. Earth- 
like planets would be easier to find around a Cen 
A than a Cen B. That is because the HZ around 
star A is more distant from this star. Naturally, 
the opposite is true for RV detections. Due to the 
difference in the stellar masses, a Cen B offers a 
better chance of finding a terrestrial planet there 
using RV techniques. The recent discovery of an 
Earth-sized planet around this star supports our 
results. The observed planetary RV signal was 
reproduced excellently by our analytic estimates 
for circular planetary orbits. 

Our prediction of RV amplitudes for terres- 
trial planets around a Cen B are also in good 



agreem ent with those presented by iGuedes et al 
(|2008l ). The four terrestrial planets used in the 



RV model by these authors produce almost exactly 
four time s the predicted RMS amplitude given in 
figure [31 iGuedes et al] (|2008l ) claim that Earth- 
like planets in the a Centauri are detectable even 
for signal-to-noise ratios (SNR) of single obser- 
vations below 0.1. However, obtaining sufficient 
data to reconstruct the planetary signal requires a 
great amount of dedicated observing time (approx- 
imately 5 years in their example). Validati ng this 
statement, it took lDumusaue et al. 1 (l2012l) about 
4 years of acquired data to detect a Cen Bb . The 
data published by iDumusque et al. ( 2012 ) also 



= 90° ±t 



t/2 llBorucki fc Summers|[l984l '). 



allow a glimpse on the current performance of the 
HARPS spectrograph revealing a precision around 
50 — 80 cm/s. Given the fact that the RV signal of 
a habitable planet around a Cen B would be still 
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half an order of magnitude smaller (Fig.|3]), consid- 
erably more observation time would be required to 
identify habitable companions. HIRES measure- 
ments are currently yielding precisions around 1 
m/s. Identiiying RV signals of habitable worlds 
around a Cen B therefore seems even more un- 
Hkely when using HIRES. The previous examples 
show that some development of observational ca- 
pacities is still necessary to achieve the RV resolu- 
tion required for discovering habitable planets in 
the a Centauri system. 

The success of NASA's Kepler space telescope 
in identifying countless Ear th-sized planetary can- 
didates (e.g.. lBoruckil 2011 ) that require follow-up 
observations might provide the necessary momen- 
tum to develop instruments capable of resolving 
RV signals in the range of cm/s. Focusing on less 
massive binaries would have the advantage of hav- 
ing greatly enhanced RV signals as the HZs will be 
situated closer to the planet's host stars. How far 
this might simplify the task of finding habitable 
worlds will be the topic of further investigations. 

In regard to transit photometry, both Kepler 
and CoRoT telescopes have proven that it is possi- 
ble to find terrestrial planets around Sun- lik e stars 



e.g.. iLeger et al.ll2009t iBorucki et al.ll2012[ ). The 



combination of proven technology and the pre- 
sented argument that the dynamical environment 
in binary star systems will enhance transit prob- 
abilities makes photometry currently the most 
promising method for finding Earth-like planets 
in the HZs of S-Type binary star systems. 

7. Summary 

In this work, we provided an analytic frame- 
work to estimate the detectability of a terres- 
trial planet using radial velocity (RV), astrometry 
(AM), as well as transit photometry (TP) in copla- 
nar S-Type binary configurations. We have shown 
that the gravitational interactions between the 
stars of a binary and a terrestrial planet can im- 
prove the chances for the planet's detection. The 
induced changes in the planet's eccentricity en- 
hance not only RV and AM peak amplitudes, but 
also the probability to witness a planetary tran- 
sit. Next to the presented "best case" estimates, 
we offered RMS/averaged expressions which are 
deemed to be more suited to determine the long- 
term influence of the second star on planetary fin- 



gerprints in S-Type systems. In contrast to peak 
amplitudes, the RMS of a planet's AM signal is 
only modified slightly by the additional gravita- 
tional interaction with the second star. A similar 
behavior can be seen in planetary transit probabil- 
ities. The RMS values of RV signals are altogether 
independent of the secondary's gravitational influ- 
ence, assuming that the system is nearly coplanar. 

After defining the Permanent, Extended, and 
Average Habitable Zones for both stellar compo- 
nents of the a Centauri system, we investigated 
the possible interaction between the newly discov- 
ered a Cen Bb and additional terrestrial compan- 
ions in a Cen B's HZ. Our results suggest that 
a Cen Bb is on an orbit with very low eccentric- 
ity which would not be influenced signiflcantly by 
habitable, terrestrial companions. Conversely, a 
Cen Bb's presence would also not affect Earth-like 
planets in the habitable zone of a Cen B. 

We estimated the maximum and RMS values of 
the RV as well as AM signal for a terrestrial planet 
in the a Centauri habitable zones. The peak and 
RMS amplitudes of the RV signal ranged between 
4 and 12 cm/s. Astrometric signals were estimated 
to lie between 1 and 5 /ias. Given the current ob- 
servational facilities, enormous amounts of observ- 
ing time would be required to achieve such preci- 
sions. If the a Centauri was a transiting system, 
however, a habitable planet could be detectable 
using current technologies. It seems that the de- 
tection of Earth-like planets in circumstellar hab- 
itable zones of binaries with Sun-like components 
via astrometry and radial velocity is still some- 
what beyond our grasp, leaving photometry to be 
the only current option in this respect. 

SE and EP-L acknowledge support from FWF 
through projects AS11608-N16 (EP-L and SE), 
P20216-N16 (SE and EP-L), and P22603-N16 
(EP-L). SE acknowledges support from the Uni- 
versity of Vienna's Forschungsstipendium 2012. 
NH acknowledges support from the NASA Astro- 
biology Institute under Cooperative Agreement 
NNA09DA77A at the Institute for Astronomy, 
University of Hawaii, and NASA EXOB grant 
NNX09AN05G. SE and EP-L would also like to 
thank the Institute for Astronomy and NASA As- 
trobiology Institute at the University of Hawaii- 
Manoa for their kind hospitality during the course 
of this project. The authors are thankful to Niko- 



11 



laos Gcorgakarakos for his valuable suggestions 
and to the anonymous referee for his/her con- 
structive comments. 

REFERENCES 

Beauge, C, Ferraz-Mello, S., & Michtchenko, 
T. A. 2007, Planetary Masses and Orbital Pa- 
rameters from Radial Velocity Measurements, 
ed. R. Dvorak, 1 

Blaes, O., Lee, M. H., & Socrates, A. 2002, ApJ, 
578, 775 

Borucki, W. J., & Summers, A. L. 1984, Icarus, 
58, 121 

Borucki, W. J., et al. 2012, ApJ, 745, 120 

Borucki, W. J., et al. 2011, ApJ, 736, 19 

Borkovits, T., Erdi, B., Forgacs-Dajka, E., & 
Kovacs, T. 2003, A&A, 398, 1091 

Chauvin, G., Beust, H., Lagrange, A.-M., & 
Eggenberger, A. 2011, A&A, 528, A8 

Correia, A. C. M., Laskar, J., Farago, F., & Boue, 
G. 2011, Celestial Mechanics and Dynamical 
Astronomy, 111, 105 

Doyle, L. R., et al. 2011, Science, 333, 1602 

Dumusque, X. et al., 2012, Nature, 491, 207 

Eggl, S., Pilat-Lohinger, E., Gcorgakarakos, N., 
Gyergyovits, M., & Funk, B. 2012, ApJ, 752, 
74 

Einstein, A., Infold, L., & Hoffmann, B. 1938, An- 
nals of Mathematics, 39, 65 

Fabrycky, D., & Tremaine, S. 2007, ApJ, 669, 1298 

Ford, E. B., Quinn, S. N., & Veras, D. 2008, ApJ, 

678, 1407 

Ford, E. B., ct al. 2012, ApJ, 750, 113 

Forgan, D. 2012, MNRAS, 422, 1241 

Georgakarakos, N. 2002, MNRAS, 337, 559 

— . 2003, MNRAS, 345, 340 

— . 2005, MNRAS, 362, 748 

Gilliland, R. L., et al. 2011, ApJS, 197, 6 



Guedcs. J. M., Rivera, E. J., Davis, E., Laughlin, 
G., Quintana, E. V., & Fischer, D. A. 2008, 
ApJ, 679, 1582 

Haghighipour, N., 2010, Planets in Binary Star 
Systems, Springer, New York 

Hestroffer, D., Dell'Oro, A., Cellino, A., & Tanga, 
P. 2010, in: Lecture Notes in Physics, ed. 
J. Souchay & R. Dvorak, Vol. 790, pp. 251, 
Springer Verlag, Berlin 

Holman, M. J., & Wiegert, P. A. 1999, AJ, 117, 
621 

Kane, S. R., & von Braun, K. 2008, ApJ, 689, 492 

Kane, S. R., Horner, J., & von Braun, K. 2012, 
ApJ, 757, article id. 105 

Kasting, J. F., Whitmirc, D. P., & Reynolds, R. T. 

1993, Icarus, 101, 108 

Kervella, P., Thevenin, F., Segransan, D., 
Berthomieu, G., Lopez, B., Morel, P., & 
Provost, J. 2003, A&A, 404, 1087 

Kiseleva-Eggleton, L., & Eggleton, P. P. 2001, in: 
Evolution of Binary and Multiple Star Systems, 
ed. P. Podsiadlowski, S. Rappaport, A. R. King, 
F. D'Antona, & L. Burderi , Astronomical So- 
ciety of the Pacific Conference Series, 229, pp. 
91 

Lee, M. H., & Peale, S. J. 2003, ApJ, 592, 1201 

Leger, A., et al. 2009, A&A, 506, 287 

Marchal, C. 1990, The three-body problem, Else- 
vier 

Ohta, Y., Taruya, A., & Suto, Y. 2005, ApJ, 622, 
1118 

Orosz, J. A., et al. 2012a, A&A, 758, article id.87 

Orosz, J. A., et al. 2012b, Science, 337, 1511 

Paardekooper, S.-J., & Leinhardt, Z. M. 2010, 
MNRAS, 403, L64 

Pilat-Lohinger, E., & Dvorak, R. 2002, CeMDA, 
82, 143 

Pourbaix, D. 2002, A&A, 385, 686 
Pourbaix, D., et al. 2002, A&A, 386, 280 



12 



Quirrcnbach, A., ct al. 2011, European Physical 
Journal Web of Conferences, 16, 7005 

Rabl, G., & Dvorak, R. 1988, A&A, 191, 385 

Ricker, G. R., Latham, D. W., Vanderspek, R. K., 
et al. 2010, Bulletin of the American Astronom- 
ical Society, 42, #450.06 

Roell, T., Neuhauser, R., Seifahrt, A., & Mu- 
grauer, M. 2012, A&A, 542, A92 

Schneider, J., Dedieu, C., Le Sidaner, P., Savalle, 
R., & Zolotukhin, 1. 2011, A&A, 532, A79 

Thebault, P., Marzari, F., & Scholl, H. 2009, MN- 

RAS, 393, L21 

Torres, G., Andersen, J., & Gimenez, A. 2010, 
A&A Rev., 18, 67 

Underwood, D. R., Jones, B. W., & Sleep, P. N. 
2003, Inter. J. Astrobio., 2, 289 

Welsh, W. F., et al. 2012, Nature, 481, 475 

Williams, D. M., & Pollard, D. 2002, Inter. J. As- 
trobio, 1, 61 



This 2-column preprint was prepared with the AAS lAlJ^X 
macros v5.2. 



13 



A. Equation of the Center 



The equation of the center providing a direct relation between the true anomaly / and the mean anomaly 
M is presented up to the 6*'' order in eccentricity e: 



g3 '5g^ ^ 

f +\2e + 1 sinAf 

■' ^ 4 96 



5e2 lie* 17e^ 



24 192 



sin[2Af] 



(Al) 



sm[3A./]+^-^)sm[4M] 



12 64 

1097 



1223 



e' sin[5M] + -^e^ sin[6Af] + 0{e') 



B. Averaging of 

The averaging integrations over M and w in equations (1171) and (jlSp were carried out as in the following; 
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The integration over M is trivial. Using the partial integration technique to integrate over w, we obtain 

27r 



1 /" 5 5 

^ / -(e2)MC0s(2w)dw = - 





(e2)M,c^ cos(2w)|f + 2(e2)M.^ / sin(2w)dtj 



27r 







Here we have used the fact that (e^ 



2ir 

M,w = (l/27r) / (e^) jvfrfw does no longer depend on w. From the definition 



of averaging given by equation ^ , we have 



fia 
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A similar procedure has been applied to derive equation 
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Fig. 1. — Top: The radial velocity signal caused by an Earth-like planet orbiting a Sun-like star with 
different eccentricities (Beauge ct al. 2007). The planet is at 1 AU with oj — 45° when e ^ 0. Bottom: The 
amplitude variations of the primary's radial velocity signal due to an Earth-like planet that is subject to 
the gravitational perturbations of a second star. Both stars are Sun-like with a separation of 20 AU and 
an orbital eccentricity of 0.5. The planet's initial orbit was circular with a semimajor axis of 1 AU. Our 
analytically estimated maximum amplitude V^""^^ is also shown. 
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Fig. 2. — Top: The maximum astrometric amplitude, — l^ixj + yj) ' \f=T,.ui=o, due to an Earth-like 

planet orbiting its Sun-like host star. The planet's orbital elements are a — 1 AU, e — g™"^^ — 0.5, w — 0, 
f2 = 111°. As shown here, the maximum distance from the origin of the coordinate system is independent 
of the system's inclination with respect to the plane of the sky (i). Bottom: Evolution of the astrometric 
signal {xf,yf) caused by an Earth- like planet in a binary star system. The planet is orbiting a Cen B at a 
distance of 1 AU. The evolution of the astrometric signal is shown for 3750 periods of a Centauri AB. Since 
the system is coplanar, the changes in orientation and shape of the projected ellipse are due to variations in 
the planet's eccentricity (e) and argument of pericenter (oj). 
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Fig. 3. — Habitability maps showing tlie maximum amplitudes and RMS of tlie RV signal of the planet- 
hosting stars for the a Centauri system. The quantity Op is the semimajor axis of the terrestrial planet and 
Ch is the binary eccentricity. The color blue shows the PHZ, the EHZ is green, yellow indicates the AHZ, and 
red means that the planet is outside of any defined HZ. The purple area denotes dynamical instability. The 
horizontal line in each panel denotes the actual eccentricity of the a Centauri binary. As shown here, 1/'"°^ 
reacts strongly to enhanced binary eccentricities (top row, curved, solid lines) whereas in contrast, {{Vr))M,uj 
is independent of the binary's eccentricity {bottom row, straight, vertical lines). The straight, vertical lines 
in the top row correspond to RV amplitudes for circular orbits {V^'^^'^). See section [5^ for details. 
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Fig. 4. — Habitability maps showing the maximum amphtudes and RMS of the astrometric signals for the 
a Centauri system. The color coding is similar to figure [S] The vertical dashed lines in the top panels 
represent regions with similar values of p'^^^'^. The curved lines in these panels show regions with similar 
pfnax amplitudes. In the bottom panels, the vertical lines represent areas of equal RMS amplitudes, ((p)) m,ui- 
One can see that planetary orbits with dynamically enhanced eccentricities can have smaller semimajor axes 
and still produce similarly high astrometric amplitudes as circular orbits which are more distant from the 
host-star. 
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Fig. 5. — Graphs of the maximum (e™""^) and RMS {{e^ m ^) values of the planetary eccentricity for different 
values of the planet semimajor axis (op) and the binary eccentricity (e;,) for a Cen A iltft) and a Cen B 
(right). The purple region denotes orbital instability. The horizontal line indicates the actual eccentricity of 
the a Centauri binary. 
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Fig. 6. — An additional terrestrial planet in a Cen B's HZ affects the orbit of a Cen Bb {top row) and 
vice versa {bottom row). In the top left panel, the numerically computed evolution of the eccentricity of a 
Cen Bb in the Newtonian three body problem (3BP) consiting of the binary a Cen AB and the planet a 
Cen Bb (blue) is compared with different four body problem scenarios (4BP). In one scenario, an additional 
Earth-sized body orbits a Cen B at the inner edge of its AHZ (red, see Table [3]). In the other scenario the 
terrestrial planet is assumed to be at the outer edge of a Cen B's AHZ (orange). The analytic estimate for 
the maximum eccentricity (e™"^) in the 3BP is presented as well (dotted, horizontal line). The top right 
panel shows the exact same setup, only with general relativity (GR) taken into account. The orbit of a 
Cen Bb becomes practically circular. While the influence of an additional planet at the outer edge of the 
HZ is barely noticable in the eccentricity evolution of a Cen Bb (blue, orange), a planet at the inner edge 
of the AHZ would cause distinct features (red). In contrast, neither GR nor a Cen Bb will influence the 
eccentricity evolution of planets (cp) in the HZ significantly as is shown in the bottom left panel. Also, the 
semimajor axes evolution of additional planets in the HZ is negligible (bottom right panel). Here, Aop/ap 
denotes the normalized difference between the 4BP+GR and the 3BP semimajor axis evolution of planets 
at the inner (red) and outer (green) edges of a Cen B's AHZ respectively. 
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Fig. 7. — Graphs of the ratio of transit probabihties (pt/ptc) with pT = PT|e=e™»=^ and pxc = PT\e=Oi in a- 
binary similar to the a Centauri system. The graphs show the transit probabihties in terms of the planet's 
argument of pericenter (w), as measured from the line of sight for component A (top) and component B 
(bottom). The red, full line corresponds to the planet starting at the inner border of AHZ and dashed- 
green corresponds to the outer edge, for each star (see Table [3]). The full and dashed straight lines are 
the corresponding ratios of the averaged transit probabilities {{pt)/{ptc)) evaluated using equation (PT|) . 
Compared to the transit probability ratios px/pxc, the ratio of averages (pt) / {ptc) shows only a weak 
dependence on the planet's initial position in the HZ. 
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Table 1 



Physical and orbital parameters of the a Centaur i ABb system (IKervella et al 
GuEDES ET AL.llioOsI : IPourbaix ET AL.ll2002l : IDumusoue et al. II2012I ). 



2003: 



a Centauri 


A B 


spectral classification 


G2V KIV 


mass \Mq\ 


1.105 ± 0.007 0.934 ± 0.007 


Teff^ [K] 


5790 5260 


luminosity [Lq] 


1.519 0.500 


distance [pc] 


1.339 ± 0.002 


period (Pb) [d] 


29187 ± 4 


ab [AU] 


23.4 ± 0.03 


eb 


0.5179 ± 0.00076 


ib [deg] 


79.205 ± 0.0041 


uJb [deg] 


231.65 ± 0.076 


VLb [deg] 


204.85 ± 0.084 


a Centauri 


B b 


Pv [d] 


3.2357 ±0.0008 




(fixed) 


minimum mass (m™*") [MEarth\ 


1.13 ±0.09 



Table 2 

A comparison between the predicted RV values using the analytic expressions derived in 

sectio n m AND THE OBSERVED VA LUES FOR THE TERRESTRIAL PLANET DISCOVERED AROUND a CeN B 



([DUMUSQUE ET AL. II2012I ). COPLANARITY OF THE SYSTEM WAS ASSUMED. ThE FORMAL 



UNCERTAINTIES HAVE BEEN DERIVED ASSUMING GAUSSIAN ERROR PROPAGATION OF THE UNCERTAINTIES 

GIVEN IN Tab. [TJ The maximum predicted planetary eccentricity for a Cen Bb is 
^max _ Q 003. When taking general relativity into account, however, the orbit of a Cen 
Bb will remain practically circular (see Fig. [6]). 





Predicted signal [m/s] 


Observed signal [m/s] 


{{yr))M,^ 


0.365 ±0.029 




TTcirc 
* r 


0.517 ±0.041 


0.51 ± 0.04 


V J. 


0.519 ±0.041 
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Table 3 

Detectability of an Earth-like planet in the HZs of the a Centauri system. Each row 

SHOWS THE maximum AMPLITUDE OF THE RADIAL VELOCITY SIGNAL AS WELL AS THE ASTROMETRIC 
FINGERPRINTS OF A TERRESTRIAL PLANET IN THE a CeNTAURI HZS. ThE CRITICAL PLANETARY 
SEMIMAJOR AXIS (Cc) INDICATES THE ONSET OF DYNAMICAL INSTABILITY ( HOLMAN fc WiEGERtI 119991 ) . 

C omputations using chaos indicat ors are in good agreement with those stability limits 
(IPilat-Lohinger fc DvoRA"i3l2002h . Analytic expressions for calcula ting the boundary 

VALUES OF PLANETARY SEMIMAJOR AXES IN THE SYSTEM'S HZS ARE GIVEN IN (IEGGL ET AL.I|2012| ). 



OL Ccn 


[AU] 


inner AHZ 


inner EHZ 


inner PHZ 


outer PHZ 


outer EHZ 


outer AHZ 










1.03 


1.07 


1.12 


1.81 


1.94 


2.06 


HZ border 


[AU] 


A 


2.76 


8.97 
5.89 


8.83 
5.78 


8.66 
5.65 


7.14 
4.44 


6.97 
4.30 


6.82 
4.17 




[cm/s] 






2.28 
1.53 


2.37 
1.59 


2.49 
1.66 


4.20 
2.69 


4.52 
2.88 


4.84 
3.06 










0.62 


0.64 


0.65 


1.13 


1.19 


1.23 


HZ border 


lAUl 


B 


2.51 


12.21 
8.25 


12.09 
8.16 


11.94 
8.05 


9.37 
6.12 


9.19 
5.98 


9.04 
5.86 


ymax 

((Vr))M,^ 


[cm/s] 






1.58 
1.09 


1.62 
1.11 


1.66 
1.14 


2.97 
1.98 


3.12 
2.08 


3.26 
2.16 


{{p))m,^ 
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